Correlation between electrons and vortices in quantum dots 
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Exact many-body wave functions for quantum dots containing up to four interacting electrons are 
computed and we investigated the distribution of the wave function nodes, also called vortices. For 
this purpose, we evaluate the reduced wave function by fixing the positions of all but one electron 
and determine the locations of its zeros. We find that the zeros are strongly correlated with respect 
to each other and with respect to the position of the electrons and formulate rules describing their 
distribution. No multiple zeros are found, i. e. vortices with vorticity larger than one. Our exact 
calculations are compared to results extracted from the recently proposed rotating electron molecule 
(REM) wave functions. 

PACS numbers: 73.21. La, 71.10.-w 



I. INTRODUCTION 

The discovery of the fractional quantum Hall effect 
(FQHE) indicated the existence of a new state of matter 
corresponding to a novel type of strongly correlated quan- 
tum many-body state. Already the first steps towards the 
understanding of this effect involved the addition of extra 
zeros to the many-particle wave function in order to ac- 
count for the electron-electron correlation. The Laughlin 
wave function^ at filling factor v = l/(2p + 1) reads 



lie* 
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where units are used such that the magnetic length is set 
equal to unity. Here, z = x — iy is a complex number 
expressing the two-dimensional coordinates of the elec- 
trons. If one fixes the coordinates of all electrons except 
one, the resulting wave function will have zeros of order 
2p + 1 located at the positions of all fixed electrons. The 
wave function Q embodies the strong correlation be- 
tween the electrons as the wave function (and the prob- 
ability to find an electron) in the vicinity of one of the 
fixed electrons vanishes more quickly than prescribed by 
the Pauli exclusion principle alone. We also note that in 
Laughlin's wave function all the zeros are rigidly bound 
to the electrons and there are no free zeros. In this re- 
spect, the wave function is rather special since a dif- 
ferent distribution of zeros (e. g., around or between the 
fixed electrons) would also be able to serve the purpose 
of stronger correlation and reduced interaction energy. 

In the subsequently formulated composite fermion 
(CF) theorji2i2£ the strong correlations were dealt with in 
a different way, by introducing weakly interacting quasi- 
particles. Also here, the zeros of the many-body wave 
function played a central role. A zero in the wave func- 
tion can also be interpreted as a vortex, i.e., when going 
around a zero its phase changes by 27m, and the winding 
number n equals the order of the zero. The new quasipar- 
ticles of the CF theory were interpreted as electrons with 



an even number of vortices or magnetic field fluxes at- 
tached to them. When a particle moves around a closed 
loop it encircles the usual Aharonov-Bohm flux due to 
the external magnetic field which is partly cancelled by 
the vortices attached to the electrons. Therefore, the 
quasiparticles can be regarded as moving in an effective 
magnetic field which is much weaker than the applied 
magnetic field. 



When constructing the CF wave function a Jastrow 
factor [zk — zi) 2p was introduced for each pair of elec- 
tron coordinates, quite similarly to Laughlin's wave func- 
tion. Subsequently, the lowest-Landau-level (LLL) pro- 
jection procedure, was introduced with the consequence 
that the vortices are no longer rigidly bound to the elec- 
trons. Thus, the relative distribution of zeros and elec- 
trons becomes less restrictive, and their correlations in 
this composite fermion liquid were investigated numeri- 
cally in recent papers 



In the present paper, we investigate the electron- vortex 
correlations in a finite system by starting from exact 
many-body wave functions obtained by means of a direct 
numerical diagonalization. Our model system is a quan- 
tum dot containing a few (up to four) electrons. The 
numerical results are compared to those obtained from 
the analytically available rotating-electron-molecul o 9 i 10 
(REM) wave functions. This recently formulated theory 
appeared as a competitor— (or at least an alternative) to 
the CF approach. It is derived from a more solid theoret- 
ical background, and introduces no a priori requirements 
on the positions of the zeros of the wave function. 



The paper is organized as follows. The model, numeri- 
cal procedure and the REM wave functions are described 
in Sec. II. The simple case of a two-electron quantum dot 
is described in Sec. III. Three- and four-electron dots are 
the subject of Sees. IV and V, respectively, and our con- 
clusions are given in Sec. VI. 
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II. COMPUTATIONAL PROCEDURE 

Let us consider a parabolic quantum dot with N = 
2 ... 4 electrons placed into a perpendicular magnetic 
field of strength B. We work with dimensionless oscil- 
lator unitsjii that is, lengths are measured in oscillator 
lengths ao = (h/m*uJo) 1 ^ 2 and energies in Huq where u>q 
is the confinement frequency. Then the relative strength 
of the inter-electron interaction is given by the dimen- 
sionless coupling constant A = ao/a* B expressed as a ra- 
tio of the oscillator length to the effective Bohr radius 
a* B = eh 2 /e 2 m*. Here, e is the dielectric constant of the 
medium and m* is the effective electron mass. The mag- 
netic field strength is expressed as a ratio of the cyclotron 
and confinement frequencies 7 = uj c /ujq. The resulting 
dimensionless Hamiltonian 



H = 



N 



i>3 



(2) 



is solved by direct diagonalization in the subspaces of 
given angular momentum L, which is a good quantum 
number. The results regarding the dependence of the 
ground state angular momentum on the confinement and 
the magnetic field strength for the case of four electrons 
in the dot were analyzed in Ref. 0. Throughout this 
paper we take A — 2 which is a typical value for exper- 
imental realized quantum dots>i£ We found that higher 
values of A require longer calculations timesA^ and did 
not result in new physics. 

In the present work we concentrate on the fully po- 
larized ground states and investigate the information 
encoded in the corresponding ground state many-body 
wave function 



*(ri 



(3) 



or, to be more precise, on the reduced wave function, 
which depends only on the position of one electron while 
the coordinates of the remaining N — 1 electrons are set 
to fixed values. Due to the Pauli principle, the reduced 
wave function has zeros at the positions of the fixed elec- 
trons. There are additional zeros which are not fixed at 
the electrons whose distribution will be the main object 
of interest in the present work. In order to locate the 
positions of the vortices we first locate the positions of 
the minima of the squared absolute value of the reduced 
wave function using a standard procedure of steepest de- 
scent from a randomly chosen initial point. Then, by 
performing a walk along a small circle around the lo- 
cated points and inspecting the change of the phase of 
the wave function we are able to distinguish actual zeros 
from other minima and determine their order, i. e. the 
winding number. 

We complement the results obtained from the exact 
diagonalization (ED) by those given by the REM wave 



functions AiS These functions are available analytically 
and help to make some exact statements. These func- 
tions are constructed by placing Gaussians at the classi- 
cal positions of electrons in strong magnetic fields and a 
subsequent restoration of symmetry. For a small number 
of electrons (N < 5), the electrons crystallize into a sin- 
gle ringi£ and the resulting wave function of the angular 
momentum L reads 10 
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Here, Zj denote the complex electron coordinates mea- 
sured in units l c y2 with l c — ^JhcjeB being the mag- 
netic length, and T> is the Slater determinant 



V(l u l 2 



,l N )= det^ 1 ,^ 2 



N 



(5) 



The wave function describes spin-polarized states of an- 
gular momentum L = Lq + Nm where L$ — N(N — l)/2 
is the smallest possible angular momentum of N spin- 
polarized electrons in the lowest Landau level, and m is 
a non-negative integer. 

Already from the general form of the REM wave func- 
tion (0J several conclusions regarding the distribution of 
zeros of the reduced wave function can be drawn. First, 
as far as the positions of zeros are concerned, the expo- 
nential factor in Eq. (JIJ can be ignored, that is zeros can 
be found from the linear combination of the Slater deter- 
minants which expands into a homogeneous polynomial 9 
Pl[z] of order L. Therefore, scaling the coordinates of all 
fixed electrons Zj — > azj ,j = 1, . . . , N — 1 results in scal- 
ing of the positions of zeros. Moreover, the polynomial 
Pl[.z] is translationally invariant 9 , therefore rigid shifting 
of the positions of the fixed electrons by the same amount 
results in a rigid translation of the distribution of zeros 
in the plane. Due to the circular symmetry of the 
quantum dot, the distribution of wave function zeros is 
also invariant with respect to rotation of the system as a 
whole. 

Regarded as a function of zn the order of the polyno- 
mial P L [z] is q N = L - (N - 1)(N - 2)12, This follows 
from the fact that in order for one electron to occupy 
the orbital with the largest possible angular momentum 
I, the remaining N — 1 electrons must reside in the low- 
est possible momentum states I = 0, 1, . . . , iV — 2. Thus, 
according to a fundamental theorem of algebra the total 
number of zeros is L for two electrons, L — 1 for three 
electrons, and L — 3 for a four electron system. 

The question of the number of zeros obtainable from 
exact diagonalization is more subtle. If the ED procedure 
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FIG. 1: The phase of the reduced wave function for two 
electrons. One electron is fixed at r = (1,0). The angular 
momentum of the state is L = 7 and we find a seventh or- 
der zero located exactly at the fixed electron. Lengths are 
measured in units of ao 



includes only the single-electron states from the lowest 
Landau level, the resulting wave function is a polyno- 
mial times a Gaussian, i. e. has a similar structure to 
the REM wave function but the expansion coefficients 
are now determined numerically. Thus, in the case of the 
lowest-Landau-level approximation we expect to find the 
same number of vortices as predicted from the REM wave 
function. However, if higher Landau levels are included 
the ED wave function (with the exponential removed) 
involves Laguerre polynomials of the argument | Zj | 2 and 
thus becomes non-analytical. This fact prevents us from 
making any exact statements regarding the total number 
of zeros. However, in the high magnetic field limit the 
lowest Landau level approximation becomes rather accu- 
rate and inclusion of higher Landau levels modifies the 
calculated wave function only at large distances from the 
quantum dot center. Thus one may still expect to find 
the same number of zeros as predicted from the earlier 
argument. 

On the other hand, the non-analyticity of the ED wave 
function makes it possible to observe besides vortices also 
anti-vortices, i. e. zeros around which the phase winds in 
the opposite direction. 



III. TWO-ELECTRON QUANTUM DOT 

For the sake of completeness, we begin with the sim- 
plest case of two electrons in a dot. We evaluate the re- 
duced wave function in a ground state of angular momen- 
tum L = 7 and plot its phase as a function of the coordi- 



nates in Fig.^ One electron is fixed at (x, y) = (1, 0)ao- 
Different shades of gray correspond to different values of 
the phase between — ir and it. Zeros of the wave function 
are located at the points where the phase is not deter- 
mined and the winding number indicates the order of the 
zero. We see that in the present case we have a single 
zero of seventh order located at the position of the fixed 
electron. 

This result can be easily understood by recalling that 
in a parabolic confinement potential the center-of-mass 
(CM) motion and the relative motion can be separated. 
The CM motion is not affected by the electron-electron 
interaction. In the ground state the CM motion is in its 
lowest state and its wave function is just a Gaussian of 
the CM coordinate which does not lead to the appearance 
of any zeros. The wave function of the relative motion at 
small values of the relative coordinate r = ri — Y2 behaves 
as ~ r m e 1 " 1 ^ where the relative angular momentum m 
coincides with the total momentum L. Therefore, in the 
reduced two electron wave function one always finds just 
a single "giant" vortex of vorticity L. 

Note that this situation is special to the parabolic con- 
finement case, and deviations from perfect parabolicity 
leads to splitting of the multiple vortex to a system of 
several single vortices. This was found in the reduced 
wave function of two electrons in a confined trion»=& where 
the non-parabolicity of the potential felt by electrons was 
due to the presence of a nearby hole. 



IV. THREE ELECTRONS IN A DOT 

The locations of the zeros calculated with the ED 
method for the N = 3 case are shown in Fig. [21 for the 
spin-polarized ground states up to L = 21. The angular 
momenta of these states are multiples of 3 as predicted 
by the magic number theory^ The two fixed electrons 
are located at (x, y) = (±l,0)a o . We observe that all 
zeros appear on a straight line defined by the two pinned 
electrons (crosses in Fig. |2J. This result persists also 
when the two pinned electrons are located off the y = 
axis, and the distances between the zeros are not visibly 
changed. 

In the case of L = 3 we find only two vortices located at 
the positions of the fixed electrons. Increasing the angu- 
lar momentum to L = 6 results in the addition of three 
vortices, one is placed between the fixed electrons and 
one on each side. The total number of observed vortices 
in this case is L — 1 as predicted by the simple estimate. 
Proceeding to higher angular momenta, we see that at 
each step one more vortex is inserted between the fixed 
electrons. Whether or not each time one extra vortex 
is added on each outer side of the electrons is difficult 
to say. The reason is that at large distances from the 
origin (typically r > 4), the accuracy of the wave func- 
tion becomes insufficient due to the limited basis set used 
in the numerical calculation. Inaccuracies may result in 
"ghost" vortices. Therefore, the calculations were limited 
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FIG. 2: The location of the zeros of the reduced wave function 
in a three-electron quantum dot for different values of the 
angular momentum L, calculated with the ED method. Two 
electrons are fixed at r = (1,0) and r = ( — 1,0). The zeros 
located on the pinned electrons are indicated by crosses and 
the dots mark the free zeros. Note that all zeros fall on a 
single straight line. Lengths are measured in units of ao. 
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FIG. 3: The same as in Fig. H but now for the REM wave 
function. The number of observed zeros is the same, however, 
their clustering around the fixed electrons is not seen. The 
arrows indicate vortices which are outside the plotted region. 
They are only indicated on the right side, but evidently are 
also present on the left side. Note that for the REM result 
lengths are originally measured in units of l c V2, but can be 
perfectly rescaled to units of ao to match with the ED result. 



to r w 3.5 (beyond which the electron density becomes 
very small, i. e. typically |^| 2 < 10~ 6 ) and some vor- 
tices located outside this region may be overlooked. On 
the other hand, for the REM method, all vortices can be 
found, including those outside r < 4, which are indicated 
with arrows in Fig. 

Addition of the new vortices between or on the outer 
side of the fixed electrons leads to a rearrangement of 
the vortices which were already present. The vortices 
are pushed towards each other and in particular to- 
wards the fixed electrons. This can be nicely seen for 
L = 9, 12, 15, 18, 21 in Fig. [3 where the pinned electrons 
(indicated by the crosses) are approached by two vortices, 
and later by four. 

With increasing 7 higher angular momentum states be- 
come the ground state. Keeping the angular momentum 
fixed and letting 7 change shows that the effect of the 
magnetic field on the positions of the vortices is surpris- 
ingly small. For example, for the L = 12 state the posi- 
tion of the vortex around 0.75 ranged from 0.74 to 0.76 
for < 7 < 15. This implies that the position of the 
vortices is mainly determined by the value of the angular 
momentum. Since we are not interested in the exact po- 
sition of the vortices but rather in their general behavior 



and their interactions, we will from now take 7 = 2. 

In Fig. 01 we show the distribution of the zeros for the 
same angular momenta and the same fixed electron posi- 
tions as in Fig. |2 but now obtained from the REM wave 
functions. We see that most of the qualitative features 
are well reproduced except for their positions in case of 
L > 6. Note that between the fixed electrons the zeros 
are more uniformly distributed and the clustering around 
the electrons as seen for the ED method is absent. 

As a matter of fact, for the case of three electrons the 
positions of zeros in the REM theory can be calculated 
exactly. Introducing the center-of-mass coordinate z = 
(zi + z 2 + 23)/3 and two relative (Jacobiii) coordinates 



z\ + z 2 



Z3 



Zb 



Z\ - Z 2 

s/2 ' 



(6) 



and dropping the exponential factors in Eq. 10} the poly- 
nomial part of the REM wave function can be written^ 
in a particularly simple form 



Ph{z a , Zb) = (z a + iz b ) L - {z a - izb) 1 



(7) 
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Equating (JJJ to zero and taking the L-th order root we 
find 

z a + izj) — (z a — izb) exp(27rifc/L), k = 1, 2, . . . , L — 1. 

(8) 

Note that there are L — 1 roots as the meaningless root 
k = has to be omitted. Eq. JSJ) is readily solved with 
the result z a — Zb cotan(7rfc/L), and using the specific 
values 2:1,2 = (±1, 0) we find the positions of the roots 



> 



cotan 



irk 

T 



k = 1,2,. ..,L- 1. (9) 



Note that despite using the specific values for the coor- 
dinates of the fixed electrons this result is still general 
since employing the above discussed symmetry proper- 
ties of the REM wave functions any randomly fixed two 
electron positions (for TV = 3) can be mapped on (±1,0). 

The result given in Eq. |J3} correctly predicts the ap- 
pearance of all zeros on a single straight line and reveals 
a simple rule for their distribution. The angular momen- 
tum for which the REM function is valid must be a mul- 
tiple of 3, that is L — 3n with n being an integer. There- 
fore, among the roots @ there always are two (namely, 
k = n and k — 2n) whose positions Z3 = ±1 coincide 
with the fixed electrons. Moreover, these two solutions, 
k = n and k — 2n, divide the interval k = 1, . . . , L — 1 
into three equal parts. Thus, we can confirm the rule 
which was already apparent in the ED results: each time 
when the angular momentum is increased by 3, three 
new vortices enter the quantum dot, and one of them is 
placed between the fixed electrons while the other two 
on the outer sides (some of the latter zeros are indicated 
outside the plotted x-region in Fig.^J). This is in agree- 
ment with our ED results when we limit ourselves to the 
\x\ < 3.5 region. Our results are in variance with those 
of Saarikoski et al. [8j who found only the addition of a 
single vortex when L increases to its subsequent allowed 
value. 

Note that the analytic expression 10 obtained from 
the REM theory fails to predict the clustering of zeros 
around the fixed electrons. Namely, Eq. @ suggests that 
the density of vortices is largest around Z3 = and mono- 
tonically decreases to both sides, this is opposite to the 
ED result which clearly shows that the density of vortices 
tends to increase around the fixed electrons and is some- 
what lower right in the middle between the two electrons. 
The REM approach is unable to reflect the subtle inter- 
action between the electrons and the zeros but at larger 
distances from the pinned electrons it predicts the zeros 
at approximately the correct positions, as can be seen for 
L = 6, 12, 18 around x = 3. 



V. FOUR ELECTRONS 

In the N = 4 case the three pinned electrons can be 
placed in many different ways. We consider three main 
configurations: a half-square triangle (corresponding to 
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FIG. 4: The location of the zeros of the reduced wave function 
for four electrons for different angular momenta L, calculated 
with the ED method. Three electrons are fixed at r = (±1, 0) 
and r = (0, 1) forming a half-square configuration. Lengths 
are measured in units of ao. 



the classical positions in a Wigner crystal), an equilat- 
eral triangle and a line configuration. 

Fig. 21 shows the positions of zeros corresponding to 
the half-square triangle configuration, calculated with 
the ED method. The pinned electrons are located at 
r = (±1,0) and r = (0,1), i.e. at the three corners 
of a square, and the considered angular momenta are 
L — 6, 10, 14, 18, i. e. the ones corresponding to the full 
spin polarization in the ground state. One immediately 
notices that the positions of the zeros of the wave func- 
tion are arranged on rays (shown by the thin lines in 
Fig. 0] Again, it is possible to spot a simple rule analo- 
gous to the one obtained for the preceding three-electron 
case that explains the location of the zeros. At the lowest 
possible angular momentum L = 6 there are three zeros 
whose positions coincide with the pinned electrons. Each 
time the ground state angular momentum is increased by 
four, four new zeros are added. One is placed inside the 
triangle defined by the three pinned electrons and the 
other three end up on the rays outside the triangle. In 
this case we looked at points up to r « 4 from the ori- 
gin, thus some of the zeros are located outside this region, 
where l^l 2 is negligible small. One notices again that the 
free zeros, seem to gather around the pinned electrons, 
which is clearly seen, e. g., for L = 18. Note that the 
L = 18 state corresponds to the v = 1/3 Laughlin state 
following the formula v = jv ^~ 1 - ) . In Fig. 0] one can 
actually see three vortices (one attached to the pinned 
electron and two free vortices) in the close neighborhood 
of the pinned electrons. 

The number of zeros inside the triangle formed by the 
pinned electrons increases by one each time the angu- 
lar momentum increases by four. So, for L — 14 there 
are two zeros located inside the triangle, and it is inter- 
esting to see how their arrangement can agree with the 
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FIG. 5: The location of the zeros of the reduced wave function 
for four electrons when L = 14 for the equilateral triangle con- 
figuration, calculated with the ED method. Three electrons 
are fixed at r = (±1,0) and r = (0, \/3). An anti-vortex ap- 
pears around r = (0, 0.5) indicated by the open dot. The inset 
shows a contourplot of the phase near the position (indicated 
by the square in the main figure) of the anti-vortex. Lengths 
are measured in units of ao. 



external symmetry defined by the pinned electrons when 
the half-square triangle is transformed into an equilateral 
one. As can be seen from Fig. [S] instead of two zeros in- 
side the triangle four zeros are formed. One of them is 
placed into the center and actually is an anti- vortex (see 
inset of Fig. [5] for a contourplot of the phase of the wave 
function), while the other three vortices are arranged on 
the vertices of an equilateral triangle, thus the symmetry 
adapts to the external symmetry and the total vorticity 
is preserved. Apparently this configuration is preferred 
over merging of two zeros into a single one with vorticity 
two (which would be sufficient to adapt to the symme- 
try). This phenomenon shows that the zeros do not like 
to sit on the same spot, and there is a certain repulsion 
between them, i. e. there is a clear tendency not to form 
vortices with winding number n > 1. As a rule we may 
state that we can expect the formation of anti- vortices by 
this rule whenever the symmetry implied by the pinned 
electrons forces vortices to come close. This will be the 
case for L = 26 in this system, but it will also be the case 
in systems with more electrons. 

This result is in contrast with the REM result which 
predicts that the two zeros inside the triangle join into 
one giant vortex, as can been seen from Fig.[S]for L = 14. 
Apparently, the REM is not capable of describing the 
subtle interaction between the zeros due to its restriction 
to analytic wave functions. One also notes that again in 
REM there is no congregation of zeros around the pinned 
electrons in Fig. H3 and like for N — 3 the vortices rather 
tend to accumulate in the center between the electrons. 

One may wonder how strong the exact location of the 
fixed electrons influences the positions of the vortices. 
Therefore, we consider the case of a line arrangement 




FIG. 6: The location of the zeros of the reduced wave function 
for four electrons for different angular momenta L, calculated 
with the REM method. Three electrons are fixed at r = 
(±1,0) and r = (0, V3). The zeros located on the pinned 
electrons are indicated with a cross and the free zeros are 
indicated with a dot. Note that the zero in the middle for 
L = 14 is in fact a giant vortex of vorticity 2 as is apparent 
from the inset which shows a contourplot of the central region. 
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FIG. 7: The location of the zeros of the reduced wave func- 
tion for four electrons for different angular momenta L, cal- 
culated with the REM method. Three electrons are fixed at 
r = (±1,0) and r = (0,0). The zeros located on the pinned 
electrons are indicated with a cross and the free zeros are 
indicated by a dot. 



of the electrons. In Fig. [7| we show the location of the 
zeros for a line configuration calculated with the REM 
method. The three electrons are fixed at r = (±1, 0) and 
r = (0,0). Notice that for L = 14 there are two vortices 
between the electrons in contrast to the previous cases 
shown in Figs. ^ and [5] where only one vortex is situated 
in the area between the electrons. The reason is that a 
single vortex would be located on top of the middle elec- 
tron. The system tries to prevent to have higher-order 
zeros and resolves this issue by taking one of the vortices, 
which was previously (see Figs. 0] and EJ) outside the in- 
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TABLE I: The polynomials Ql(z) for the case when three 
pinned electrons are situated on a single line. 



L 


Ql(z) 


6 


1 




10 


3z 4 - 18z 2 + 10 




14 


32 s - 52z 6 + 212z 4 - 448z 2 


+-80 


18 


3z 12 - 102z 10 + 990z 8 - 6160z 6 4 


- 14003z 4 




-7837z 2 +291 





ner electron region, and placing it between the electrons 
which results in a symmetric distribution of vortices. In 
contrast to the two electron case, there are zeros which 
appear next to the line defined by the three pinned elec- 
trons. When we look at the locations we can derive a 
simple rule that explains the addition of the vortices: ev- 
ery time one goes to the next magic angular momentum 
four zeros are added. The first time they are added on 
the y = line and are equally distributed in between 
the pinned electrons. In the next step they are added 
symmetrically above and below the y = line. These 
two rules alternate each time the angular momentum is 
increased by four. 

In the case of four electrons in the dot it is not possi- 
ble to derive and solve a general compact expression for 
the polynomial describing the distribution of zeros in the 
REM wave function as it was done for the three-electron 
dot. For the present configuration featuring the arrange- 
ment of three pinned electrons into one line such poly- 
nomial has the form z(z 2 — I)Ql(z) where the first two 
factors represent the zeros located on the fixed electrons, 
and the polynomial Ql(z) describes the distribution of 
the vortices. In Table 1, we give this polynomial for the 
considered values L = 6, 10, 14, 18. Note that thanks to 
the symmetry of the configuration only even z powers 
appear in Qi(z). 

Next we compare the previous REM results with our 
exact calculation. Therefore, we show in Fig. [HI the same 
configuration as in Fig.0but now the results are obtained 
with the ED method. The location of the zeros is quali- 
tative similar to those for the REM functions, but again 
we see that there is a much stronger clustering around 
the fixed electrons and the vortices above and below the 
y — line are much further away. 

To summarize the dependence of the location of the 
vortices on the positions of the fixed electrons we show 
in Fig. [5| a 3D plot for L = 14 in which we fix two elec- 
trons at r = (±1,0) and move the third electron along 
the vertical axis from r = (0, 0) to r = (0, 2). Notice that 
we can clearly see: (1) how the vortices move with chang- 
ing symmetry of the fixed electron distribution, (2) the 
appearance of an anti- vortex when the two inner vortices 
come close to each other and (3) how the anti-vortex ex- 
ist over a certain range before merging with a vortex and 
being annihilated, and (4) how the position of the two 
vortices are rotated over 90° with respect to the position 
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FIG. 8: The same as in Fig. but calculated with the ED 
method. Three electrons are fixed at r = (±1,0) and r = 
(0,0). Lengths are measured in units of ao 




FIG. 9: The position of the inner vortices as a result of the 
movement of the fixed electrons for L = 14. Two electrons 
are fixed at r = (±1,0) and the third one is fixed at r = (0, z). 
The projections of the vortex positions on the x — and y = 
planes are given by dashed lines for clarity. The triangle itself 
is indicated by dotted lines. The open dots indicate the region 
of existence of the anti-vortex. Only the three fixed electrons 
forming the triangle and the vortices inside the triangle are 
shown. Lengths are measured in units of ao. 
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dence of repulsion between the zeros and their attraction 




FIG. 10: Same as in Fig.HHbut now calculated with the REM 
method. Notice that there is no anti- vortex present. 

of the electrons when one electron moves away from the 
two others, i. e. with increasing z. This rotation occurs 
through the intermediate creation of an anti-vortex. 

Similar results for the REM reduced wave function are 
shown in Fig.^3 One observes that also here the relative 
position of the vortices are rotated over 90° but that in 
this case no anti-vortex is formed to make this happen, 
i. e. the two vortices join in one giant vortex of vortic- 
ity two and then they separate again into two distinct 
vortices. 

Another interesting thing to investigate is the evolu- 
tion from a triangle configuration towards a line config- 
uration. As Fig. 03 shows for L — 14 two vortices are 
located above the y — line at r = (±2.6, 1.8) and two 
below at r = (±2.6,-1.8). When we start moving the 
fixed electron at r = (0, 0) upwards the two vortices lo- 
cated above the y = line will move closer to each other 
and finally a vortex-anti-vortex pair will be created as 
soon as they are close enough as shown in Fig. ^] Going 
further a vortex and an anti- vortex will meet and annihi- 
late, changing the configuration to one with two vortices 
located on the x = line. The mechanism behind this is 
the same as shown in Fig. [!|] 



3.5 
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FIG. 11: Vortex positions calculated with the ED method 
for L — 14. Three electrons are fixed at r = (±1,0) and 
r = (0, 0.273). The zeros located on the pinned electrons are 
indicated with a cross and the free zeros are indicated by a 
dot. A anti-vortex (open dot) appears at r = (0,3). The 
inset shows the phase of the reduced wave function around 
the anti- vortex. Lengths are measured in units of do- 



to the pinned electrons leading to a strong correlation 
between the vortices and between the electrons and the 
individual vortices. Additional vortex-anti-vortex pairs 
can be formed for certain symmetries of the fixed elec- 
tron distribution. Qualitatively, several of the results on 
the distribution of the zeros can be obtained from an an- 
alytically available rotating-electron-molecule wave func- 
tions. However, the REM theory is not able to describe 
the condensation of zeros around fixed electrons and the 
formation of an anti-vortex. 



VI. CONCLUSIONS 

We investigated the distribution of zeros of the reduced 
many-body wave function in few-electron quantum dots. 
The results show that the arrangement of zeros can be 
described by a set of simple rules. The number of vortices 
increases with Al = N between two subsequent fully po- 
larized magic angular momentum states. The vortices (or 
zeros) between the electrons are situated on rays point- 
ing away from the electron cluster. There is clear evi- 



Acknowledgments 

This work is supported by the Belgian Interuniversity 
Attraction Poles (IUAP), the Flemish Science Founda- 
tion (FWO-V1), VIS (BOF) and the Flemish Concerted 
Action (GOA) programmes. EA. was supported by the 
Marie Curie fellowship program under contract number 
HPMF-CT-2001-01195. 



9 



* Electronic address: maarten.tavernier@ua.ac.be 

t Electronic address: egidijus.anisimovas@ua.ac.be 

* Electronic address: francois.peeters@ua.ac.be 

1 D. C. Tsui, H. L. Stormer, and A. C. Gossard, Phys. Rev. 
Lett. 48, 1559 (1982). 

2 R. B. Laughlin, Phys. Rev. Lett. 50, 1395 (1983). 

3 J. K. Jain, Phys. Rev. Lett. 63, 199 (1989). 

4 J. K. Jain, Phys. Rev. B 41, 7653 (1990); ibid. 42, 9193(E) 
(1990). 

5 Composite Fermions: a unified view of the quantum Hall 
regime, editted by O. Heinonen (World Scientific, Singa- 
pore, 1998). 

6 K. L. Graham, S. S. Mandal, and J. K. Jain, Phys. Rev. 
B 67, 235302 (2003). 

7 D. Pfannkuche and A. H. MacDonald, Physica E (2004). 

8 H. Saarikoski, A. Harju, M. J. Puska and R. M. Nieminen, 
cond-mat/0402514 

9 C. Yannouleas and U. Landmann, Phys. Rev. B 66, 115315 



(2002) . 

10 C. Yannouleas and U. Landmann, Phys. Rev. B 68, 035326 

(2003) . 

11 M. B. Tavernier, E. Anisimovas, F. M. Peeters, B. Szafran, 
J. Adamowski, and S. Bednarek, Phys. Rev. B 68, 205305 
(2003). 

12 S. Tarucha, D.G. Austing, T. Honda, R. J. V. der Hage, 
and L. P. Kouwenhoven, Phys. Rev. Lett. 77, 3613 (1996). 

13 S. A. Mikhailov, Phys. Rev. B 66, 153313 (2002). 

14 W. Y. Ruan, Y. Y. Liu, C. G. Bao, and Z. Q. Zhang, Phys. 
Rev. B 51, 7942 (1994). 

15 V. M. Bedanov and F. M. Peeters, Phys. Rev. B 49, 2667 
(1994). 

16 E. Anisimovas and F. M. Peeters, Phys. Rev. B 68, 115310 
(2003). 

17 L. Jacak, P. Hawrylak, and A. Wojs, Quantum Dots 
(Springer, Berlin, 1998). 



